Institute for Space Studies 



PRICE $ 


I PRICE(S) $ 


Norman R. Simon 


ard copy (HC) 


and 


licrofiche (MF)__ Richard Stothers 

Ju ' v65 Institute for Space Studies 

Goddard Space Flight Center, NASA 


New York 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 



PULSATIONAL INSTABILITIES IN HYDROGEN- POOR 


MASSIVE BLUE STARS, I. 


Norman R, Simon 


and 


Richard Stothers 

Institute for Space Studies 
Goddard Space Flight Center, NASA 
New York 



ABSTRACT 


The structure and pulsational properties of massive stars 

with helium cores and thin, hydrogen-poor envelopes have been 

investigated. The core structures are very insensitive to 

modifications of the hydrogen envelope,. However, it is found 

that the total stellar radius and luminosity are modified by 

the integrated hydrogen content and by the hydrogen gradient, 

respectively. The models are pulsationally stable under modest 

conditions of central condensation (hydrogen-burning shell and/or 

high hydrogen content) . The unstable models are energized almost 

entirely by the helium reactions in the core, while most of the 

damping occurs in the thin, non-burning hydrogen envelope. For 

each mass, the period (~ hour) is remarkably insensitive to the 

envelope modifications. The critical models dividing stable 

from unstable models are determined by (1) a pulsational eigen- 
2 

frequency w ~ 3 and (2) the value of the hydrogen gradient. 

Our critical models explain some of the observed features of 
the classical Wolf-Rayet stars of high luminosity but seem to be 
too blue with respect to Rublev’s sequence of stars on the H-R 


diagram. 



2 


I. INTRODUCTION 

Massive stars on the main sequence become unstable against 
nuclear-energized pulsations above a certain critical mass 
(Ledoux 1941; Schwarzschild and Harm 1959). However, evolution 
off the main sequence quickly stabilizes stars up to the greatest 
masses (Schwarzschild and Harm 1959). Even when the energy source 
shifts to a shell outside the contracted helium core, blue stars 
persist in exhibiting high pulsational stability (Ledoux 1941; 

Cox 1955; Stothers and Simon 1968). The relevant physical argu- 
ments have been amply reviewed in the foregoing references. 

All calculations of models so far have been based on the 
assumption that the stars evolve at constant mass and with devel- 
oping chemical inhomogeneity. When either or both of these 
assumptions are relaxed, pulsational instability becomes a greater 
possibility. As a trivial formal example, consider the evolution 
of a completely mixed star during hydrogen burning. As hydrogen 
is depleted throughout the star, the critical mass for instability 
decreases as the inverse square of the mean molecular weight 
(Schwarzschild and Harm 1959) ; 

M crit /M 0 ■ 21 /<' 2 - 

A combination of somewhat more complex factors may result 
in more physically realistic examples. It is therefore of interest 
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to ask what observational evidence exists for stars of high 
mass and blue color (because of the necessarily small central 
condensation) showing a significant degree of instability. The 
classical Wolf-Rayet stars are certainly prime candidates. 
However, their evolutionary stage is still a mystery and has 
been the subject of much discussion. Various proposals have 
been guided by observational indications placing the stars in 
different positions with respect to the main sequence, as well 
as by the relatively low observed masses in binary systems, by 
some very high luminosities, and by the apparent chemical 
peculiarities of the stars. 

Among the more noteworthy proposals are the following. Pre- 
main sequence : gravitationally contracting stars (Sahade I960; 

Underhill 1966). Main sequence : stars at the end of hydrogen 

burning (Westerlund 1961, 1964) ; stars having completely mixed 
at the end of hydrogen burning and now undergoing substantial 
mass loss (Westerlund and Smith 1964) ; pulsationally unstable 
stars in secondary contraction at the end of hydrogen burning, 
in analogy with the 8 Cephei stars (Stothers 1965). Post-main 
sequence : stars evolved beyond hydrogen burning with substantial 

loss of mass (Limber 1964; Rublev 1965; Tanaka 1966; Paczynski 
1967); stars having completely mixed before, during , or after 
helium burning (cf. Burbidge, Burbidge, Fowler, and Hoyle 1957); 
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and stars in the helium-burning stage (Salpeter 1953; Crawford 
1953; Divine 1965) . Clearly, virtually every possible evolution- 
ary stage has been suggested! 

In the present paper, we propose to study the structure 
and pulsational properties of post-main-sequence stars with 
massive helium cores and thin, hydrogen-poor envelopes. Pulsa- 
tional instability (energized by nuclear reactions) may account 
for the observed variations in many Wolf-Rayet stars. The very 
high luminosities of some of the stars are assumed in the present 
work to be the consequence of high mass. Nevertheless, the 
original masses must have been still greater in order for the 
stars to have attained the assumed configurations. 


f 
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II. ASSUMPTIONS 

The stellar models under consideration contain three zones: 
(1) a radiative outer zone with a smooth gradient of hydrogen 
and helium extending from the surface down to a thin shell 
(where hydrogen may or may not be burning) ; (2) a radiative 

intermediate zone, devoid of hydrogen, below the shell; and 
(3) a convective core where helium burning occurs. 

The hydrogen abundance in the radiative outer zone is 
assumed to be a linear function of mass fraction 

X = a Q + a 1 q (1) 

where a Q and a^ are determined by fixing the hydrogen content 
X_ at the surface (q = 1) and requiring X = 0.03 at the shell 

S 

(q ) . The shell is approximated by a discontinuous jump in 
s 

hydrogen and helium as well as in luminosity. For simplicity, 
we shall consider only the initial model of helium burning in 
the core, so that the helium abundance below the shell may be 
represented by Y = 1 - Z, where Z is the metals abundance. 
Deinzer and Salpeter (1964) and Boury and Ledoux (1965) have 
shown in the case of stars composed of pure helium that little 
change in the structure and pulsational properties takes place 
as helium is depleted in the core. On this basis we assume that 
our results would not be altered much by the choice of a more 
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evolved core. The metals abundance is taken to be Z = 0.03 
throughout our models. 

The assumptions adopted here are similar to the ones 
adopted by Stothers (1966a) and by Stothers and Simon (1968) , 
except that the present models lack a zone of zero-age chemical 
composition capping the zone of changing composition. The 
opacity is assumed to be due solely to electron scattering: 


x = 0.19 (1 + X) (2) 

The nuclear-energy release is provided by the CNO-cycle in the 
shell and by the triple-alpha reaction in the core. The rate of 
the former is given by 

e == e° XX e T V H (3) 

H H CNO r 17 


with 


X CN0 “ Z//2 , V H : 


14, log e H = -99.0, 


while the rate of the latter is given by 


He 


0 3 2 V H <a 


(4) 


The quantities and are allowed to vary from model to 

model since they are strongly dependent on temperature (Reeves 


1965) . 
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III. STABLE MODELS 

The masses selected for calculation are 15, 60 and 100 M® . 

The other free parameters which have to be specified are X R and 

q . The method of calculation has been described in Stothers 
s 

(1966ai) . Results for each mass using a variety of envelope 
parameters are given in Table 1. 

For a given mass, it is found that the luminosity remains 

approximately constant while the spatial distribution of matter 

in the star (indicated by the stellar radius or by the radial 

extension of the shell) varies with the positioning of the 

shell in mass fraction q and with the surface abundance of 

s 

hydrogen X . For increasingly homogeneous models (q “* 1), 

Jtx S 

both the total luminosity and the fraction of luminosity con- 
tributed by hydrogen burning decrease. Furthermore, the radius 

fraction of the shell increases with q and the total radius 

^•s 

decreases. Similar results were found by Cox and Salpeter 
(1961) for models of stars with much lower masses. 

It is useful to differentiate two ways of formally altering 
the hydrogen content in the envelope: (1) addition of hydrogen^ 

holding the shell mass fraction fixed and (2) reduction of the 
shell mass fraction / holding the surface hydrogen abundance 
fixed. In both cases, the integrated hydrogen content J X dq 
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is increased and so is the stellar radius . However, the 
hydrogen gradient dX/dq is steepened in case (1) and reduced 
in case (2) ; the stellar luminosity (as well as the hydrogen- 
burning fraction of the luminosity) decreases and increases, 
respectively. 

For models with a fixed core mass, alterations of the 
envelope composition have virtually no effect on the structure 
of the core. Thus, the core is effectively "decoupled" from 
the envelope and behaves like an isolated star (cf. Hayashi, 

Hoshi, and Sugimoto 1962) , 

A comparison of the various masses having the same envelope 

parameters (X , q ) shows that the fraction of luminosity con- 

tributed by hydrogen burning is smaller at the higher masses. 

Depending on the selected envelope parameters (particularly q ) 

s 

this fraction varies from less than 1 percent to nearly 40 percent 
in our models. However, for any constant set of envelope para- 
meters, the basic stellar quantities show the same trends with 
mass as do normal main-sequence stars. 

On account of the small amount of mass contained in the 
hydrogen envelopes, our models lie considerably to the left of 
the core hydrogen-burning main sequence (Stothers 1966b) , 
although, of course, not as far left as the analogous main 
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sequence for pure helium stars (Deinzer and Salpeter 1964) , 
Increasing the mass of the hydrogen envelopes still further 
can result in location of the models in the blue-supergiant 
or even the red-supergiant region (Stothers and Chin 1968; 
Stothers and Simon 1968) . 
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IV. PULSATING MODELS 

The radial pulsation characteristics of our models were 
calculated in the usual linear, quasi-adiabatic approximation 
(Schwarzschild and Harm 1959; Boury and Ledoux 1965; Stothers 
and Simon 1968) . The relative pulsation amplitudes are all 
continuous through the star, except for the luminosity perturba- 
tion which has a discontinuity at the hydrogen shell. The 
average rate at which energy is being fed into, or removed 
from, the pulsations over a period is given by 

P PN PH PS (5) 

where I. is the rate of gain from nuclear sources, L is the rate 
PN W 

of flux damping, and L_ c is the rate of damping by acoustical 
waves running off the surface. A positive sign for L indicates 
an unstable configuration. The energizing term L pN contains a 
contribution both from shell hydrogen burning L pN and core helium 
burning L ptJ , while the main damping term L pH is the sum of the 
heat loss rate in the envelope L® and in the convective core 

Ext 

L p . To avoid the difficulty involved in calculating luminosity am- 
plitudes in a convective region (cf. Boury, Gabriel, and Ledoux 
1964) , we have estimated them following Schwarzschild and Harm 
(1959) by means of an average involving quantities computed at 
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the center and outer boundary of the core only. Such an approxi- 
mation is sufficient for our purpose because the contribution 
from these terns is small, never exceeding 30 per cent of the 
contribution from the surface waves. 


The calculated pulsational rates are listed in Table 2 along 


with the square of the dimensionless frequency 


U) 


( — ) 
V Period / 


_R_ 

GM 


( 6 ) 


Period 

the pulsational _e - folding time 1/K, where K is the stability 
coefficient (L p / 2E p ) and E p is the mechanical energy of the 
pulsations,* and the relative pulsation amplitudes at the surface, 
shell, and center. The quantities with asterisks have been com- 
puted without the contribution from surface waves. Otherwise, the . 
notation is the same as in previous papers ( Schwarz schi Id and 
H&rm 1959; Stothers and Simon 1968). 


The general conclusions that can be drawn regarding the pul- 

2 

sational eigenfrequencies are already known. For example, u) 

increases with greater central condensation (as measured by p /p ) 

c 

and decreases with greater relative radiation pressure (as measured 

by 1 - b ) • For reference, we note (1) a greater central conden- 
c 

sation with increasing radius and (2) a greater relative radiation 
pressure with increasing luminosity and hence with increasing mass. 
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The critical value of iw dividing the stable and unstable models 

2 

is approximately the same for all masses (to ~ 3 ) . The reason 

is that the critical central condensation increases along with 

relative radiation pressure at higher masses (see below and 

2 

Appendix) . This result for the critical values of uu seems to 
be valid not only for our hybrid models but also for pure helium 
stars (Boury and Ledoux 1965) and for Population I hydrogen 
main-sequence stars ( Schwarz schi Id and Harm 1959). 

The pulsation periods of those four models lying in the 

vicinity of the critical model for each mass are remarkably 

similar. This result is surprising in view of the fact that the 

integrated hydrogen abundance and the hydrogen gradient (with 

respect to mass fraction) vary by factors up to 6 and 8, respect- 

2 

ively. The increase of uu with radius is apparently such that 
the period is held nearly constant (eq. [6] ). 

m all the models, most of the damping comes from heat 
leakage in the envelope. The core contributes less than 5 per 
cent of the interior damping and the surface waves contribute 
less than 20 per cent of the total damping (interior plus surface) . 
Energizing of the pulsations comes almost entirely from the 
helium reactions in the core. The hydrogen shell contributes 
less than 10 per cent of the energy for pulsations in all the 
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unstable models. This is in general due to the fact that L /L 

H 

is small in the unstable models, although the smallness of L /L 

H 

is not a sufficient condition for instability. When the shell 
is vigorously burning and/or hydrogen is added to the envelope, 
the envelope is in an expanded state. Such an expansion drives 
up the central condensation because the central density is 
fixed at a nearly constant value by the helium-burning reactions. 
Consequently, the pulsation amplitudes drop off more rapidly 
inside the surface (see Appendix) and therefore are smaller near 
the center where most of the destabilization takes place. Hence, 
a high central condensation tends to induce stability. The 
critical value of p c /p dividing the stable and unstable models for 
each mass is given in Table 3. The corresponding value of 6r/r 
at the center is about 40 per cent of the surface amplitude in 
all cases. 

In Table 3 we have also listed the critical value of q for 

^s 

each selected mass and surface abundance of hydrogen. An asterisk 
again indicates that surface waves have been neglected. In 
general, the deeper the mass fraction of the shelj|., the less 
surface abundance of hydrogen is required to induce stability. 

The key factor determining the critical model for a given mass 
seems to be the hydrogen gradient rather than the total amount 
of hydrogen. This is so since the hydrogen gradient determines 
the degree of central condensation. On account of the smaller 
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radiation pressure at lower mass, the low-mass models become 
unstable under more restrictive conditions on the hydrogen 
gradient. In the limit of vanishing hydrogen content (pure 
helium stars) the critical mass is 7 to 8 M 0 (Boury and 


Ledoux 1965) . 
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V. COMPARISON WITH OBSERVATIONS 

The location of Wolf-Rayet stars on the H-R diagram is 
very uncertain and there is much confusion and argument concern- 
ing where they lie with respect to the main sequence. This is 
true as regards both luminosity and effective temperature. 

In addition, the masses of these stars are not well known. 

There are, however, some generally agreed upon properties of 
Wolf-Rayet stars which are of interest here, although they will 
be discussed in more detail in a subsequent paper. 

First of all, the classical Wolf-Rayet stars are rare objects 
and must be young because of their close association with 0 and 
B stars. Second, these objects seem to be super luminous (as 
compared with main-sequence stars) for any reasonable masses that 
are assigned them. Third, many of the stars are surrounded by 
expanding shells and often by dense, spherically symmetric nebulae. 
Fourth, their emission lines are frequently observed to vary 
irregularly on a time scale of hours. Fifth, their electron 
temperatures are very high and their effective temperatures may 
also be very hot (possibly hotter than main-sequence stars) . Sixth, 
the atmospheres of these stars seem to be rich in helium and, in 
the case of the WN sequence, possibly rich in nitrogen as well 
(as if the surface material had at one time been processed through 
the GN-cycle) . 
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It has often been suggested in the literature that Wolf- 
Rayet stars may arise as a result of mass loss. Rublev (1965) 
states: "It cannot be excluded that WR stars are some sort of 

remnants of initially very massive stars that have shed a con- 
siderable portion of their matter (including all of the 
hydrogen-rich envelope) and have uncovered layers which duririg 
the earlier stages of evolution corresponded to the periphery 
of the convective core," 

The foregoing evidence allows us to make the tentative 
suggestion that some Wolf-Rayet stars may be massive objects 
in the core helium-burning stage of evolution and, further, 
that the structure of these stars could be approximated by our 
hydrogen-poor , helium-burning models . 

Let us examine this suggestion. Figure 1 shows a plot of 
the H-R diagram of our model stars (open circles), including 
pulsationally unstable models (crossed open circles) . The 
filled circles represent Wolf-Rayet stars as located by Rublev 
(1965). The pure helium sequence (Deinzer and Salpeter 1964) 
and the normal main sequence (St others 1966b) along with approx- 
imate mass locations are drawn in for reference, at left and 
right respectively. 
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It is immediately clear that all Wolf-Rayet stars on the 
diagram are superluminous , the least luminous star corresponding 
to a main-sequence mass of about ,60 M . This suggests that the 
Wolf-Rayet stars are not hydrogen-burning, but represent rather 
a post-ma in- sequence phase of evolution. If such an evolutionary 
stage is indeed reached via mass loss on the upper main sequence, 
then the expanding shells and nebulae which surround Wolf-Rayet 
stars may be attributed to the ejected matter. 

As previously noted, Wolf-Rayet masses are extremely uncertain. 
Underhill (1966) estimates the average mass in binary systems as 
lying between 4 and 8 M while., if one takes Rublev's stars to 
be burning core helium, his highest mass appears to be ~ 90 M 0 , 

Let us turn now to the question of Wolf-Rayet stability. 

As can be seen in Figure 1, a narrow band of instability exists 
for our hydrogen-poor models, running close to the pure helium 
sequence and turning off slightly at high masses. Since Rublev's 
estimates of effective temperature are higher than most quoted 
in the literature, it seems unlikely (though still possible) 
that the Wolf-Rayet locus can be moved any farther to the left. 

In that case helium core-energized pulsations cannot be the 
cause of observed variability, since the cooler of our models 
prove to be quite stable. 
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It is clear that any mechanism which will account for the 
present observed properties of luminous Wolf-Rayet stars should 
also fit in with the previous evolutionary history of mass loss. 
In a forthcoming paper, we shall propose just such a unified 
mechanism. 


One of us (N.R.S.) gratefully acknowledges the support of 
a Yeshiva University fellowship during the academic year 1966- 


67 . 



APPENDIX 


It is required to show that the gradient of the radius 
amplitude at the surface is steeper for a higher central 
condensation. The relevant linearized equations for adiabatic 
pulsations are 



At the surface, x = 1, q = 1, and V = - dlnP/dln r = 00 . Therefore, 
we need the regularity condition 





The gradient at the surface is clearly 



setting (6r/r) = 

R 




1 . 



We now make use of the well-known variational result (Ledoux 


and Walraven 1958) 

uu 2 £ (3? - 4) j 

where the ratio 



is a measure of the central condensation, like g /p . Clearly, 

o 

2 . 

ou increases with increasing central condensation and with 

decreasing relative radiation pressure. As an approximation, we 

set r = T and use the equality sign in the variational expression 
R 

for uo 2 . Then 



At the stellar surface is determined from 


S R = 1 “ (H r L/ 4ttcGM) . 


For a given mass and chemical composition, the surface value 
of g and thus of T depends only on the luminosity. But the 
luminosity is not very sensitive to chemical composition changes 



in the envelope. Hence, for a given mass. 


6r 
r 

indicating the more rapid drop-off of pulsation amplitudes 
(all of which depend on 6r/r) with higher central condensation. 
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TABLE 1 


Model Characteristics for Helium-Burning Stars 
with Hydrogen-Poor Envelopes. 


M = 15 M@ 


X R 

0.10 

0.10 

0.05 

0.05 

g s 

0.95 

0.98 

0.95 

0.97 

8 

s 

0.681 

0.708 

0.680 

0,712 

log T 

s 

7.651 

7.518 

7.654 

7.577 

log c s 

0.876 

0.532 

0.882 

0.715 

CO 

0.563 

0.661 

0.590 

0.655 

q 4 

0.763 

0,786 

0.763 

0.785 


0.592 

0.591 

0,592 

0.591 

log T 
■ c 

8.287 

8.287 

8.287 

8.287 

log p c 

2.637 

2.636 

2.637 

2.636 

p c/p 

34.4 

27.7 

29.8 

25.0 

V L 

0.089 

0.000 

0.110 

0.004 

log (L/L ) 

© 

5.471 

5.440 

5.481 

5.441 

log (R/R ) 
© 

0.075 

0.044 

0.054 

0.029 

log T 

e 

5.093 

5 . 101 

5 . 106 

5.109 

v He 

20.0 

20.0 

20.0 

20.0 

i 0 

log e He 

-164.82 

-164.82 

-164.82 

-164.82 







TABLE 1 (continued) 


M = 60 M* 
© 


0.05 


0.15 

0.80 

0.200 

7.751 

0.245 

0.125 

0.664 

0.372 

8.327 

2.368 

3560. 

0.377 

6.411 

1.037 

4.847 

18.0 

-148.23 


0.15 

0.95 

0.372 

7.685 

0.421 

0.487 

0.879 

0.347 

8.333 

2.339 

66.8 

0.045 

6.359 

0.471 

5.118 

17.0 

-139.90 


0.15 
0.98 
0.388 
7 . 557 
0.067 
0.589 
0.922 
0.347 
8.334 
2 . 340 
50.8 
0.000 
6.355 
0.431 
5.136 
17.0 
-139.90 


0.10 

0.95 

0.368 

7.689 

0.426 

0.524 

0.876 

0.347 

8.333 

2,339 

53.6 

0.069 

6.370 

0.439 

5.136 

17.0 

-139.90 


0.90 

0.247 

7.743 

0.338 

0.332 

0.770 

0.354 

8.331 

2.346 

201 . 

0.299 

6.438 

0.628 

5.059 

18.0 

-148.23 


0.05 

0.95 

0.374 

7.695 

0.453 

0.570 

0.879 

0.346 

8.333 

2.338 

40.6 

0.075 

6.374 

0.399 

5.157 

17.0 

-139.90 




table 1 (continued) 


M = 100 M© 


0.15 

0.95 

0.289 

7.691 

0,274 

0.463 

0.894 

0.279 

8.348 

2.245 

83.7 

0.033 

6.636 

0.609 

5.118 

17.0 

-139.90 


0.15 

0.97 

0.298 

7.620 

0,080 

0.523 

0.924 

0.278 

8.349 

2.246 

69.2 

0.001 

6.635 

0.581 

5.131 

17.0 

-139.90 


0.10 

0.90 

0.197 

7.742 

0.209 

0.287 

0.796 

0.284 

8.346 

2.251 

329. 

0.222 

6.680 

0.805 

5.031 

17.0 

-139.90 


0.10 
0.95 
0.291 
7.698 
0.300 
0.517 
0. 895 
0.278 
8.348 
2.245 
58.9 
0.046 
6.642 
0.558 
5.145 
17.0 
-139.90 


0.05 

0.90 

0.189 

7.748 

0.206 

0.340 

0.792 

0.284 

8.346 

2.249 

196. 

0.245 

6.692 

0.731 

5.070 

17.0 

-139.90 


0.05 

0.93 

0.237 

7.732 

0.283 

0.454 

0.842 

0.280 

8.347 

2.246 

82.8 

0.165 

6.676 

0.607 

5.129 

17.0 

-139.90 







TABLE 2 


Pulsational Characteristics for Hydrogen- Poor Models. 

(Asterisks indicate neglect of running waves.) 


M = 15 H e 


X R 

— - - 1 

0.10 

0.10 

0.05 

0.05 

q s 

0.95 

0.98 

0.95 

0.97 

2 

uu • 

4.044 i 

3.403 

3.546 

3.076 

Period (hr) 

0.46 

0.45 

0.46 

0.45 


0.679 

0.701 

0.686 

0.713 

log p r 

-7 . 149 

-7.082 j 

-7.064 

-6.999 

log g R 

5.464 

5.525 

5.505 

5.555 

(6r/r) R 

1.000 

1.000 

1.000 

1.000 

(&T/T) r 

-2.205 

-2.050 

-2.075 

-1.971 

(6l/L) r 

-4.280- 

-4.199 

-4.300 

-3.885 

(6r/r) g 

0.416 

0.566 

0.495 

0.606 

(6T/T) s 

-0.630 

-0.887 

-0.739 

-0.921 

( 6L/L) 

S 

-0.510 

-0.887 

-0.607 

-0.876 

(6r/r) c 

0.320 

0.422 

0.384 

0.467 

(6T/T) 

C I 

-0.369 

-0.486 

-0.443 

-0.539 

L S /L 
PN / 

0.69 

0.00 

0.89 

0.06 

L C /L 
PUT 

3.199 

6.085 

4.504 

7.448 

L /L 

par 

3.89 

6.09 

5.40 

7.50 

l ph //l 

6.613 

5.416 

5.591 

4.843 

l ph //l 

0.121 

0.219 

0.174 

0.269 

l ph /l 

6.734 

5.636 

5.765 

5.112 

L P S / L 

1.129 

1.105 

1.092 

1.122 

Lp/L 

' -3.97 

-0.65 

-1.46 

1.27 

Lp/L 

-2.85 

0.45 

-0.37 

2.39 

1/K (yr) 

-725 

1 -7860 

-2730 

+4930 

1/K* (yr) 

-1010 

+11300 

I -10100 

+2620 


TABLE 2 ( continued ) 


M = 60 M 

0 


0.15 

0.15 

0.15 

0.10 

0.05 

0.05 

0.80 

0.95 

0.98 

0.95 

0.90 

0.95 

6.521 

3.505 

2.828 

2 . 866 

5.507 

2.235 

5.07 

0.98 

0.95 

0.97 

1.34 

0.96 

0.268 

0 . 350 

0.356 

0.363 

0.290 

0.386 

- 8.680 

- 7.700 

- 7.632 

- 7.590 

' - 7.936 

- 7.453 

4.140 

5.274 

5.353 

5.338 

4.959 

5.417 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

- 2.652 

- 1.906 

- 1.735 

- 1,746 

- 2.401 

- 1.590 

- 6 .610 

- 3.624 

- 2.940 

- 2.985 

- 5.603 

- 2.359 

0.009 

0.369 

0.524 j 

0 . 480 

0.088 

0.608 

- 0.014 

- 0.473 

- 0.682 

- 0.604 

- 0.129 

- 0.751 

- 0.001 

- 0.186 

- 0.355 

- 0.245 

- 0.026 

- 0.315 

0.005 

0.301 

0.417 

0.395 

0.060 

0.508 

- 0 . 006 

- 0.320 

- 0.443 

- 0.419 

- 0.064 

- 0.539 

0.00 

0.22 

0.00 

0.37 

0 . 09 

0.63 

0.000 

2.262 

4.487 

3.767 

0.069 

6.190 

0.00 

2.48 

4.49 

4.13 

0.16 

6.82 

10.534 

4.482 

3.374 

3.462 

8.234 

2.538 

0.000 

0.047 

0.113 

0.081 

0.001 

0.135 

10.534 

4.530 

3.488 

3.542 

8.235 

2.672 

0.271 

0.511 

0.450 

0.474 

0.533 

0.439 

- 10.80 

- 2.56 

0.55 

0.12 

- 8.61 

3.71 

- 10.53 

- 2.04 

1.00 

0.59 

- 8.08 

4.15 

- 0.858 

-488 

+4170 

+17400 

- 12,2 

+885 

- 0.880 

-610 

+2290 

+3460 

- 13.0 

+792 


TABLE 2 ( continued ) 


M = 100 M Q 


0.15 

0.15 

0.10 

0.10 

0.05 

0.05 

0.95 

0.97 

0.90 

0.95 

0.90 

0.93 

, 3.295 

2.847 

5.855 

2.401 

4.655 

3.015 

1.26 

1.23 

1.86 

1.23 

1.61 

1.30 

0.263 

0.265 

0.219 

0.285 

0.235 

0.262 

- 7.879 

- 7.835 

- 8.215 

- 7.720 

- 8.024 

- 7.785 

5.219 

5.275 

4.826 1 

5.322 

4.974 

5.223 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

- 1.839 

- 1.726 

- 2.477 

- 1.616 

- 2.178 

- 1.768 

- 3.354 

- 2.903 

- 5.908 

- 2.464 

- 4.710 

- 3.071 

0.368 

0.465 

0.054 

0.533 

0.143 

0.394 

- 0.453 

- 0.576 

- 0.077 

- 0.642 

- 0.193 

- 0.488 

- 0.140 

- 0.215 

- 0.010 

- 0.210 

- 0.035 

- 0.134 

0.307 

0.384 

0.038 

0.452 

0.106 

0.323 

- 0.321 

- 0.402 

- 0.040 

- 0.473 

- 0.111 

- 0.338 

0.15 

0.01 

0.02 

0.36 

0. 16 

0.62 

2.282 

3.696 

0.03 

4.900 

0.213 

2.231 

2.43 

3.71 

0.05 

5.27 

0.37 

2.85 

4.053 

3 . 349 i 

8.895 

2.721 

6.415 

3.609 

0.037 

0.069 

0 . 000 

0.081 

0.003 

0.033 

4.090 

3.418 

8.895 

2.802 

6.417 

3.641 

0.442 

0.403 

0.492 

0.391 

0.486 

0.431 

- 2.10 

- 0.12 

- 9.33 

2.07 

- 6.53 

- 1.22 

- 1.66 

0.29 

- 8.84 

2.46 

- 6.05 

- 0.79 

-537 

-14800 

- 5.69 

+1140 

- 32.4 

-999 

-680 

+5930 

- 6.01 

+957 

i 

- 35.0 
! 

-1540 
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figure caption 


Fig. 1. - Theoretical H-R diagram for super luminous massive stars. 

Filled circles are Wolf-Rayet stars as located by Rublev. 
Open circles represent our models, including unstable mo- 
dels (crossed open circles) . The stars are bounded on the 
left and right by the pure helium sequence and normal main 
sequence, respectively; several masses (solar units) are 


indicated. 
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